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Abstract 



Scattering of electromagnetic waves by many small particles of arbitrary shapes 
is reduced rigorously to solving linear algebraic system of equations bypassing the 
usual usage of integral equations. 

The matrix elements of this linear algebraic system have physical meaning. 
They are expressed in terms of the electric and magnetic polarizability tensors. 
Analytical formulas are given for calculation of these tensors with any desired 
accuracy for homogeneous bodies of arbitrary shapes. An idea to create a "smart" 
material by embedding many small particles in a given region is formulated. 

1 Introduction 

Wave scattering by small particles was studied by Rayleigh, starting in 1871. He under- 
stood that the main term in the field, scattered by a small particle, is the dipole radiation. 
A particle is small if A >> a, where a is the characteristic dimension of the particle. The 
particle is assumed homogeneous with parameters e, \i and a. By k = ^ we denote 
the wave number in the medium, surrounding particles. There is a large literature on 
scattering by bodies, small in comparison with the wavelength (see [I], [3], 0], [5] and 
references therein). Exact analytical solutions were found for spherical and ellipsoidal 
particles (jlj). 

For particles of arbitrary shapes the author gave analytical formulas for S-matrix, 
which allow one to calculate this matrix with any desired accuracy (0). Our aim in 

MSC: 78A45, 78A48, 81V10. PACS: 0200, 0340K, 0380. 

Key words: Electromagnetic wave scattering, small particles of arbitrary shapes, many-body scat- 
tering problem, nanotechnology. 



this paper is to show that the many-body scattering problem for N small particles in 
electromagnetic (EM) wave theory context can be rigorously reduced to solving linear 
algebraic systems (las) of equations bypassing the usual usage of integral equations. 

The matrix elements of this las have physical meaning: they are expressed in terms 
of the polarizability tensors of small bodies. The decisive point is: the author has de- 
rived analytical formulas which allow one to calculate these tensors with any desired 
accuracy for bodies of arbitrary shapes. Our theory is developed in Section |21 The basic 
assumption is: 

a « A << d, (1) 

where a is the characteristic dimension of the small particles and d is the smallest distance 
between two distinct particles. In Section |2] the problem is formulated and a method of 
its solution is developed. Our theory in Section |2] uses some ideas, similar to the ideas in 
[2|. The principal difference between the results in |2] and in Section |2] is the following: 
the scattering coefficients in [2] are not known analytically, and should be calculated 
separately, while in our theory the analogs of these coefficients, the scattering matrices 
Si, are given analytically, explicitly (see formulas (A.1)(A.5)-(A.7) in the Appendix). In 
[2| the isotropic point scatterers are considered, the scattering is assumed isotropic, and 
in Section |21 this assumption is not used. However, this difference is less important than 
the principal difference, mentioned above. 

In Section Elthe scattering in a medium consisting of many small particles is discussed 
under assumption At the end of this Section we formulate an idea of creating a 
"smart" material nanotechno logically, by embedding many small particles in a given 
region in such a way that the resulting material would, for instance, have a desired 
radiation pattern. 

In Section 0] the field near the boundary of small particles is discussed. In Appendix 
some auxiliary results from [S] are given. These results are essential for our theory. 



2 Statement of the problem 

Consider the scattering of a monochromatic plane EM wave with frequency u by N small 
homogeneous particles Di, 1 < % < N, with parameters e,fi,a (permittivity, magnetic 
permeability, conductivity) and Lipschitz boundaries Si, placed in a medium with pa- 
rameters Eo,{1q,<Jq = 0. Let S := U^Si, D := Uf =x Di. The time-dependent factor e~ lult 
is omitted. 

For simplicity let us assume that a is so large that condition [Af, E] = on 5* holds, 
where M is the exterior unit normal to 5* and [A/", E] := M x E is the cross product. 
Our theory can be developed for impedance boundary condition [A/", E] = ([[N, H],J\f] 

( >\ 1/2 

as well, where ( is the surface impedance, ( = ( — J , and e := e + ^. In this case the 

penetration depth 5 of the EM field into the particle is given by the formula: 5 = y^j) 
and 5 « a if ua is sufficiently large, \x being fixed, E ~ j-H if - >> e. 
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The governing equations are 

V x E = kjfiH, V x H = -iue'E in D, 



(2) 



V x E = iujfi H, V x H = -iue E in D' , (3) 
[J\f, E] and e'E ■ Af are continuous across S, (4) 

[Af, E'} = - [Af, E ] on S, E = E + E', H = H + H', (5) 

where E , H is the incident field, which satisfies (J2J) in the whole space, E', H' is the 
scattered field. Let 



and «Sj be the scattering matrix corresponding to a small particle Di. This means that if 
an electromagnetic wave IA is incident upon Di, then the scattered field W in the far- field 
zone is 

e ikr / y \ 

W = - — 5jW + o - , r = la: - xA » A >> a, (7) 

where Xi is a point in Di, which can be chosen arbitrarily, and it does not matter which 
point is chosen since Di is small. Formula (jJJ) is valid if the field U is practically ho- 
mogeneous within the distances of order a. The foundation of our theory is an explicit 
formula for S\ n ^ for a small homogeneous body of arbitrary shape, which allows one to 
calculate Si with any desired accuracy: 



; (n) 



S t \ = 0{q n ), < q < 1, 



where q is a constant which depends only on the geometry of Si and the material constants 
£j,/ij,(7j. This formula is derived by the author (0 Chapter 7]) and the results we use 
in this paper are presented below, in the Appendix. 

If Si, 1 < i < N are known and assumption ((TJ) holds, then the EM field at any point 
x G R 3 , such that 

min \x — xA > d, (9) 

l<i<N 

can be calculated by the formula: 

e %k\x-y\ 

U(x) = U (x) + V g(x,Xi)SiV(xi), g( x ,y) = f . (10) 



8=1 



The vectors V(xj) in ()10j) are unknown. The expression g>(x, Xj)<SiV(a:j) is the field, 
scattered by i-th particle, placed in the field V(xj). If V(xj) were known, then formula 
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(Jl(Jj) would give the solution to the iV-body EM wave scattering problem in the region 
0. 

If x is near the boundary Sj, then one gets: 



/ 



U{x) = Uo(x) + ^2g(x,Xi)SiV(xi) + 



V x / g(x,s)j{s)ds 
Js 3 

x V x / g(x,s)i(s)ds 



\ 



In (II lj) the integrals represent the electric and magnetic fields generated by the j-th 
particle in an immediate neighborhood of this body, and j is an unknown tangential field 
on Sj, representing the induced by U surface current. 

If assumption (JTJ) holds, one may consider the j-th particle as being placed in the 
homogeneous incident field 



g(xj,Xi)SiV( 



;i2) 



and the corresponding scattered field is 



e 

k\x — x 



g(x j ,x i )SiV(x i 



(13) 



Therefore we get a linear algebraic system of equations for the unknown V(x«) in (|T0|) : 
V{xj) = U (xj) + Y^9( x ii ^i)SiV{xi), 1<j<N. (14) 

If Sj, 1 < j < N, are known, then (|14|) is a linear system of 6N equations for the 
6N unknowns V{xj), since V(xj) is a 6-component vector. System (fPfj) can be solved 
efficiently by iterations, provided that its matrix is diagonally dominant. The matrix of 
system (fTlj) is 

I — Sig (xj ,Xj), (15) 

where / is the unit matrix in the 6-dimensional space of vectors with complex-valued 
coordinates and the norm of a vector is defined as follows: 



\\U{x)\\ := max \Ui(x)\, U 



Kt<6 









(u 4 


H 




|. H=[ 













(16) 



Matrix (|15|) is diagonally dominant if 



^2Sig(x j} Xi)\\ < 1. 



(17) 



Inequality (|T7|) holds if 

kd 



1 6 

— max y^y^i(5i) m ^i < 1. 

kd l<m<6 ^ ^ 
i=i i+j 



Here we have used the estimate of the norm of a matrix A, corresponding to the norm 
(HE): 



A\\ = max > \A m A 

Km<6 ^— ' 



fed - 



Using formula (A.l) (see Appendix) in (fTSj) one gets \(Si) m eg(xj, Xj)\ = 0( 
O(^) < 1 if assumption (JTJ) holds. For example, if ka < 0.1, kd = 20, N < 100, 
then the left-hand side of (JTHJ) is less than ^ 10" 3 6 • (N - 1) < 3 • 10~ 2 < 1, so that 
condition (fTHj) holds, and the system (fTlj) can be efficiently solved by iterations: 

Vj n+1) = U 0j + ^gjiS i Vf\ U 0j := U (xj), g j{ := g(xj,Xi), (19) 
V 7 (0) =W 0:; -, vf } :=V^(^)- (20) 



This completes the discussion of the scattering EM waves by N small particles. In this 
discussion the field on the surface of i-th particle was not calculated. 

This field is not needed for calculating the vectors U{xj) which define the total field by 
formula (fTU|) at all points except immediate neighborhoods of the particles. See Section 
4 for the discussion of the field in an immediate neighborhoods of the particles. 



3 Scattering in a medium consisting of small parti- 
cles 

If the number iV of particles is very large, say iV ~ 10 23 , so that one discusses the EM 
wave scattering in a medium consisting of small particles, then assumption Q implies 
that the number of small particles per unit volume of the space is O (^-) and their total 

volume per unit volume of the space is O (^^J — ^ if ^ — ^ 0. Therefore the sum in (jlOj) 

in this limit tends to zero if assumption holds. Assumption was essential to our 
arguments, so, if s, a and u are fixed, we arrive at the conclusion that scattering by iV 
small particles, under the assumptions (0) and in the limit liniTv^oo ^ = 0, is vanishing 
in this limit. However, the physical situation is dramatically different if we assume that 
e = s{uj) and a may depend in a suitable way on a and d. 

For example, assume that the small particles are identical spheres of radius a, a = 0. 
It is known that the polarizability tensor of a ball equals to 

— 66q — — Oij, 
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where ^ is the Kronecker symbol, i.e., the unit matrix. Recall that the polarizability 
tensor of a dielectric body D with the dielectric constant e, placed in a constant 
electrical field E in a homogeneous medium with the dielectric constant Eq, is defined by 
the equation: 

Pi = aijVeoEj, 

where V is the volume of D, over the repeated indices summation is understood, P is 
the induced dipole moment acquired by the body D in the field E, and Pj is the i-th 
component of vector P. 

Let us assume that the following limit exists: 

Vo\B(y,b) ' mPJ ' (21) 
Vol J B( 2 /,6) = ^vr6 3 , 

where the limit in (J21|) is taken as the number v of particles in the ball B(y, b) of radius 
b centered at a point y in D tends to infinity in such a way that assumption (JTJ holds 
and ^ — > 0, and /3 is the scattering direction, (3 = t|^|t if the scattering wave is directed 
from point y to point x. 

Equation (jlOj) in the limit 

N^oo, ~^0, a«X«d, k = -^-, (22) 
yields an integral equation for U(x): 

U(x)=U (x)+ [ f.^ y \ q(y,x)U(y)dy. (23) 
Jd k\x-y\ 

If the particles are identical spheres of radius a, we have 

v ~ ^-nb 3 /d 3 , and = 3e £ , ,5° kj, 
3 e + 2e 

so the nonzero limit (j21j) exists if the limit 



3 

lim °- :=w^0 (24) 



exists. This limit exists and w = £_£a if 



a 3 



e = -2e + k— , k ^ 0, (25) 

where k is a constant. There also exist materials with e = e(uj) < and o ^ for which 
the limit exists. For example, if e := e(u) + i^, a ^ 0, and 

icr a 3 

e(u) — > — 2e , ► «i -™ as — >■ cj , «a 7^ 0, (26) 
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where u is some frequency and K\ is a constant, then the nonzero limit (|24|) exists and, 
therefore, the nonzero limit (|21j) also exists. 

Physically, it is natural to assume that the self-consistent (effective) field in a medium, 
consisting of many (N ~ 10 23 ) small particles, should not change if any single particle 
is removed from the medium, in spite of the fact that the exact field can be changed 
very much in an immediate neighborhood of the removed particle ([7]). For instance, if 
the boundary condition on the surface Sj of this particle was [J\f, E] — then, after the 
particle is removed, [TV, E] 7^ on Sj, so the relative change is infinite. 

If the limit (j21|) exists, then equation ()23|) is the equation for the self-consistent field 
in the medium consisting of many small particles. 

Our result allows one to formulate an idea for creating nanotechnologically a "smart" 
material by embedding many small particles in a given bounded region. Such an idea was 
first proposed for acoustic wave scattering in see also [TO]. The idea can be briefly 
described as follows. If the limit (}2"Tj) exists, then the integral equation (J2"3|) holds. If 
for x in the far zone q(y,/3) '■— q(y,x) = q(y), so that q does not depend on x, then 
equation (}2l?|) is equivalent to the Schrodinger equation with potential q. The "smart" 
material, which is characterized by a desired scattering amplitude (radiation pattern), 
can be constructed by first solving the inverse problem of finding such a q from the 
desired scattering amplitude, that is, by solving an inverse scattering problem, and then 
relating this q to the density of distribution of the small particles, embedded in our 
region. If N(y) is the number of the small particles per unit volume in a neighborhood 
of the point y, then q(y) is proportional to N(y) (see |9J. Thus, if q is found, then N(y) 
is known, and the small particles should be embedded with the density N(y) in order 
that the resulting "smart" material would have the desired radiation pattern. Such is 
the situation in acoustic wave scattering by small bodies. The details are discussed in 
this case in [§| and [TUj . 

In EM wave scattering the potential q, in general, may depend on both variables, 
y and x. This leads to a new inverse problem of finding such a q = q(y, x) from the 
scattering data, an open problem currently. Therefore, it is of interest to find if there 
are physically reasonable conditions under which q(y,x) = q(y), because the inverse 
scattering problem in this case has been solved by the author (see [TT], chapter 5). 

4 Equation for the field in an immediate neighbor- 
hood of a particle 

Formula (jl(Jj) was derived for x satisfying condition Q. If one is interested in the field 
near Sj, one has to find j in equation and then use formula to calculate the 
field at a point x near Sj. 
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The electrical field near Sj is 



E(x) = Eo(x)+y2Ei(x)+Vx / g(x,s)j(s)ds, 

where the vector Ei(x) is calculated as the first three components of the vector 

SiU{Xi)g{Xj, Xi), 



(27) 



see formula (JTTJ). The boundary condition [Af, E]\ 
equation for j: 



and formula (|27|) yield an 



where 



Aj+j= "27 



41 



Af,E {s)+Y,Ei 



s e Sj, 1 < j < N, 



ik\s— 1\ 



' 2nk\s — 1\ 



,m 



dt, 



s G Sj. 



(29) 



Formula (|2*Hjl is analogous to formula (1.27) in [5], being its generalization to the poten- 
tials of single layer with the density j which is a tangential to Sj vector field. 

Equation (|2~KJ) is a Fredholm-type integral equation if Sj is a a sufficiently smooth 
surface. It allows one to find j uniquely if k 2 is not an eigenvalue of the homogeneous 
Maxwell's equation in D with the boundary condition 



[Af,E}\ 



0. 



If k 2 is such an eigenvalue, then this integral equation can be modified so that the modified 
equation is uniquely solvable for j (see |BJ). 

If j is found, then formula (fTTj) allows one to calculate U = 



of Sj and on Sj. 



in a neighborhood 



Appendix 

In this Appendix we collect some of the results from [5] used in Section EJ 

The formula for the iS^-matrix for EM wave-scattering by a single small body is (j3J 
p.114]): 



_ k 3 V f noP u + a 2 2 cos 9 - a 32 sin 9 a 2 i cos 9 - a 31 sin 9 - fMoPn 
^ L 47r \ai2 - cos 9 + /i ^3i sin 9 a u + ^0/^22 cos 9 - /io/%2 sin 9 

In (A.l) Se differs by the factor k from the 5-matrix in j5, p. 114] because we use the 
dimensionless function while in 15] the function ^— is used. The notations are the 

kr ' ! — - r 

same as in [5]: 
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9 is the angle of scattering, = £%(7 e ) is the electric polarizability tensor, 7 e := 
Aj = ce,j ( — 1 ) is the magnetic polarizability tensor, V is the volume of the small body, 
and Se is defined by the formula for the scattered electrical field E' in the far zone: 

ikr 

E' = —S E E, r:=\x\, (A.2) 
kr 

where E is the incident field at the point where the small body is located. If E' is known, 
then 

c p n n «^ 



H' = iufio— ik[x°,E'l x° = r- r (A.3) 



Thus, if one knows <Se, then one can calculate the matrix S in (J2J. 

The matrix <Sg in (A.l) is calculated in the coordinate system in which the wave 
vector k, |k| = k, of the incident wave and vector k' of the scattered wave are lying in 
one plane, 9 is the angle between k and k', the x'-axis coincides with the x-axis, y'-axis 
makes an angle 9 with y-axis , ^'-axis is directed along k', and the origin is the same 
for x,y,z and x',y',z' coordinates. The origin lies inside the small body (particle), and 
the plane yOz = y'Oz' is called the scattering plane. One has E = (Ei,E 2 ,0) in xyz 
coordinates, E' = (E[,E 2 ,0) in x'y'z' coordinates, and 

E 2 \ _ ( E' 2 



s *{e;) = {k)- (a - 4) 

Let us give formulas (see pp. 54-55]) for calculating the polarizability tensors 0^(7), 

7 := (e - e )(e + e ) := 7 e , 

of a homogeneous body D with boundary S, volume V, and permittivity e, placed in the 
medium with permittivity sq. We have 

l«!; } (7) - ay(7)l = 0(q n ), < q < 1, (A.5) 

where the n-th approximation cqj(j e ) to otij(%) is given by the formulas: 



m=0 v 7 



and 



b if : = / / dsdtNi^Njis) [ ■■■[ -^—i;(t 1 ,t)i;(t 2 ,t 1 )...i;(t m _ 1 ,t m ^)dt 1 ...dt m _ 1 , 



W,*)--= Wt -, r st :=\s-t\, 
4 1) (7 £ ) = «(7 £ + 7 £ 2 )^-J.^ ) . 
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By Afj(s) we denote the j-th component of the exterior unit normal to S at the point 

seS. 

The magnetic polarizability tensor is: 

:= Oy(-l) + Oijftn), 7 M := ^— — . (A.7) 

/i + /io 

If fi = jiQ then 0^(7^) = 0. 

The induced electrical dipole moment P on the body D, placed in a homogeneous 
electrostatic field E, is 

Pi = a^VeoEj. (A.8) 

Here and below summation is understood over the repeated indices. 
The induced magnetic moment is 

Mi = PijVuoHj. (A.9) 

The first term in (A.7) is absent if the penetration depth 5 of the magnetic field into D 
is much larger than the size a of D, and is present if 5 << a, for example, if the body D 
is a very good conductor in which H cannot penetrate. 

Finally, let us explain the remark below (A. 2) about the iS-matrix defined in (JJJ). 
Matrix Se in (A.l) has been obtained in [3] from formula (7.95) in [S]. This matrix 
allows one to calculate the scattered field E' by formula (A. 4). If E' is known, then the 
scattered field H' can be calculated from the Maxwell equation @: 

which is formula (A. 3), because k = uj^/Sqiiq. Thus, the knowledge of the matrix (A.l) 
allows one to calculate iS-matrix in (|ZJ). 

Let us discuss the role of the assumption ((H). We have 



e ik\x-y\ g |X| V M x 



I I V |x| \x\ 2 



X 



If \y\ « \x\, then 

Ak\x\ 
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e 



k\x\ \ \\x\ 

If 



k 


\y 


12 




x\ 





*(-°,») i + o M +O ^ , k\y\«l, \y\ « \x\. 



then 



2tt 

\y\ < Q>> F — d, d » a, a << A, A = — , 

k 



01M +0 M?)«i. 
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Thus, if |x| > d, \y\ < a, a << d and a << A, then 



\x\ 



The assumption d » A is not used in these estimates. Therefore, on the wavelength 
A there can be many small particles as long as each of these particles is in the far zone 
with respect to all other particles. This was used in [7] in a study of the acoustic wave 
scattering in a medium consisting of many (N ~ 10 23 ) small, acoustically soft, particles. 
In EM scattering one has to use the relation 

e ikr / 1\ x -y 

v x g = -j— [ik i r, r=\x-y\, 

kr \ r J \x — y\ 

and in the far zone one neglects the term - compared with ik. This can be done if 

k » -, r >d, 
d 

and these inequalities imply 

d » A. 

That is why we have imposed assumption when dealing with EM wave scattering. 
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